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ABSTRACT 

Close encounters between galaxies are expected to be a viable mechanism, as predicted by numerical 
simulations, by which accretion onto supermassive black holes can be initiated. To test this scenario, 
we construct a sample of 562 galaxies (M* > 2.5 x 10 10 M Q ) in kinematic pairs over the redshift 
range 0.25 < z < 1.05 that are more likely to be interacting than a well-matched control sample of 
2726 galaxies not identified as being in a pair, both from the zCOSMOS 20k spectroscopic catalog. 
Galaxies that harbor an active galactic nucleus (AGN) are identified on the basis of their X-ray 
emission (Lo.5-10 koV > 2 x 10 42 erg s _1 ) detected by Chandra. We find a higher fraction of AGN in 
galaxies in pairs relative to isolated galaxies of similar stellar mass. Our result is primarily due to an 
enhancement of AGN activity, by a factor of 1.9 (observed) and 2.6 (intrinsic), for galaxies in pairs 
of projected separation less than 75 kpc and line-of-sight velocity offset less than 500 km s _1 . This 
study demonstrates that close kinematic pairs are conducive environments for black hole growth either 
indicating a causal physical connection or an inherent relation, such as, to enhanced star formation. 
In the Appendix, we describe a method to estimate the intrinsic fractions of galaxies (either in pairs 
or the field) hosting an AGN with confidence intervals, and an excess fraction in pairs. We estimate 
that 17.8^7 4% of all moderate-luminosity AGN activity takes place within galaxies undergoing early 
stages of interaction that leaves open the question as to what physical processes are responsible for 
fueling the remaining ~ 80% that may include late-stage mergers. 

Subject headings: galaxies: active — galaxies: interactions — quasars: general — X-rays: galaxies 
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1. INTRODUCTION 

There has been a long-standing question in astro- 
physics, ever since quasars and A GN were firml y be- 
lieved to be powered by accretion ([Salpe tcr 1964]) onto 
supermassive black holes (SMBH); what physical mech- 
anism^) is (are) responsible for the loss of angular mo- 
mentum of gas that is initially rotationally-supported 
and t hat can then be available to fuel a central SMBH 
(e.g., iLvnden-Bell fc Reesl [l97lh ? There have been a 
number of proposed mechanisms that have been inves- 
tigated, over the past few decades, such as disk instabil- 
ities, galaxy mergers, or supernova-driven winds, just to 
name a few. More recently, galaxy mergers, including in- 
teractions between galaxies with larger separations, have 
been a leading contender si nce numerical simulations 
demonstrate that such events (jMihos fc Hernquistll 19961 ) 
can generate mass inflow to the nuclear region thus po- 
tentially fueling bot h AGNs and central starbursts (e.g. 
Hop kins et al 



qg bot l 
. 2008) 



In support of the merger scenario, the by-products of 
major mergers of gas-rich galaxies, the ultra-luminous in- 
frared galaxies (ULIRGs), have been known for decades 
to be conducive to black hole growth (jSanders Sz Mirabel 
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1996). ULIRGS have been shown to have high nuclear 
concentrations of gas (jScoville et al.l Il989| ) that likely 
provide an ample fuel reservoir for accretion onto a 
SMBH. Indeed, the inc idence of AGN in ULIRGs (e.g., 
iKartaltepe et al.l l2010t ) rises substantially with IR lu- 
minosity to an unprecedented level (> 50%) for the 
brightest that reach bolometric luminosities equivalent 

to quasars. 

There has been much eff ort (e.g.. iBahcall et all 119971 ; 
iCanalizo fc Stocktonl |200"1 over the past decade to de- 
termine whether galaxy interactions and mergers play 
a role in fueling quasars by carrying out imaging cam- 
paigns of their host galaxies with HST to look for signs 
of morphological disturbances. Such efforts are chal- 
lenging because of the glaring light of a nuclear point 
source, the lack of adequate control samples, and the 
fact that the detection of such signs of interac tion is 
highly redshift dependent (|Kampczvk et al.ll2007t) . For- 
tunately, new samples of AGNs have been constructed 
that appear to alleviate some of these difficulties. In 
particular, X-ra y surveys are identifying an obscured 
population (e.g. iBarger et alj |2003t ISzokolv et all 12004 
ISilverman etatl |2010) that makes it easier to discern 
the morphology of their host. Large-area surveys, with 
HST cover age and multi-wav e length support such a s 
COSMOS (|Scoville et alJl20Q7t iKoekemoer et al.l 12007ft . 
are enabling the construction of adequate control sam- 
ples. To date, studies based on X-ray selected AGN 
have demonstrated that major mergers of galaxies are 
not likely to be the single dominant mechanism respon 

i ggerin g 
iGabor et al.ll200l lC 
markable examples of a binary AGN have been found in 
the C OSMOS field (IComerford et al.ll2009t ICivano et al " 



sible for triggering AGN activity (iGrogin et al.l 120051 : 
]isternas et al.ll2011l ). even though re- 



201C ) and elsewhere (e.g. jGreen et al. 1 120101 iShen eTal 



20101 : iFu et al.ll201QT ) . One explanation may be that deep 



surveys mainly detect moderate-luminosity AGNs (i.e. 



"Seyfert" population, Lx- 



10 



42-44 



ergs s x ) thus 



the fueling m echanisms may differ from t he more lumi- 
nous quasars (jHopkins fc Hernquistl 120091 ). 

We present an alternative test of the hypothesis that 
galaxy interactions play a role in fueling moderate- 
luminosity AGNs. We identify galaxies in the zCOSMOS 
spectroscopic survey that are in close physical proximity 
(i.e., galaxies in close kinematic pairs). These are two 
galaxies within a given projected separation and line- 
of-sight velocity difference that likely indicates an inter- 
acting system that may merge at subsequent times. This 
approach is independent of galaxy morphology thus unbi- 
ased in this respect. We determine the fraction of galax- 
ies in kinematic pairs that host an AGN based on X-ray 
emission detected by Chandra and compare to a well- 
defined control sample. It is important to recognize that 
the X-rays enable us to cleanly identify black hole accre- 
tion, less hampered by star formation that can severely 
impact optical emission-line studies, especially for galax- 
ies in pairs. We highlight that we are investigating the 
impact of interactions on scales between ~ 10 — 150 kpc 
thus our study will provide constraints on early-stages 
of merger-driven models of black hole growth before a 
coalescence phase. Conclusive evidence for such an en- 
hancement of moderate-luminosity AGN activity, rela- 
tive to the non-active galaxy population, h as been elusive 
in SDSS samples at lower redshifts (e.g.. iLi et al.l [20081 : 



Ellis on et aill2008ft . Throughout this work, we assume 
H = 70 km s" 1 Mpc" 1 , il A = 0.75, n M = 0.25, and AB 
magnitudes. 

2. ZCOSMOS 20K GALAXY REDSHIFT SURVEY 

zCOSMOS dLillv et al.ll2007l 120091) is an ESO program 
to acquire spectroscopic redsh ifts for a large sam ple of 
galaxies in the COSMOS field (|Scoville et alj|2007t) with 
VIMOS on the VLT. A 'bright' sample of 20,682 galaxies 
(i < 22.5) is observed with a red grism to provide a wave- 
length coverage of 5500 — 9500 A ideal for identifying L* 
galaxies up to z ~ 1.2. We use this catalog to construct 
a well-defined sample of galaxies. Specifically, we iden- 
tify 15,807 galaxies with iacs < 22.5 and spectroscopic 
redshifts between 0.25 < z < 1.05, each having a quality 
flag (> 1-5) that amounts to a confidence of ~ 99% in 
the redshift measurements for the overall sample. We do 
not include any serendipitous objects that happen to fall 
within a slit. Full details on data acquisition, reduction, 
and r e dshift measurements can be found in iLillv et al.1 
(|2007l f2009h . 

We further isolate 10,964 zCOSMOS galaxies that fall 
within the Chan dra survey area (see below for details, 
lElvis et "al"1l2009ft since we use the X-rays from this data 
set to identify those hosting an AGN. Of these, there 
are 3481 with a stellar mass above M* > 2.5 x 10 10 M Q 
and 0.25 < z < 1.05. These limits are chosen to ensure a 
fairly complete representation of both bl ue and red galax- 
ies ove r the given range in redshift (see ISilverman et al.l 
2009b) . Derived properties such as stellar mass and rest- 
frame color ( Figure [J) are de t ermin ed f or each galaxy 
as det ailed in iBolzonella et all ([2010ft and lPozzetti et all 

<f2?TToh . 

2.1. Kinematic pairs 

A major focus of the zCOSMOS survey is the study 
of galaxies in close kinematic pairs as a laboratory to 
deter mine the redshift ey olution of the galaxy merger 
rate (|de Ravel et al.l 120111) and the physical prop erties 
of the galaxies themselves (jKampczvk et al.l 1201 lft . We 
use the zCOSMOS 20k catalog to identify those galax- 
ies with a nearby neighbor. Galaxies in kinematic pairs 
are specifically selected to have a projected separation 
dr < 100 kpc h _1 (143 kpc for our chosen h = 0.7) 
and a line-of-sight velo city difference dV < 500 km s _1 
(|Kampczvk et al.l 1201 lh . Based on repeat observations 
of 569 zCOSMOS galaxies, individual redshifts are accu- 
rate to a z = 0.00036 x (1 + z) or a z = 110 km s" 1 (see 
ILillv et al.ll2009l for details), an accuracy well below our 
chosen velocity difference for identifying pairs. Based on 
this selection, we identify 753 galaxies, out of the 3481 
(selected by stellar mass and redshift) that are associ- 
ated with a galaxy pair. Of these, there are 330 with a 
mass ratio less than 3 to 1, while 562 have a mass ratio 
between 10 and 1. In Figure [TJ we show the stellar mass 
and rest-frame color of the galaxies in kinematic pairs. 

3. CHANDRA X-RAY OBSERVATIONS AND AGN 
IDENTIFICATION 

We match our zCOSMOS galax y sample to the 
catalog of Chandra X-ray sources ([Elvis et~aT1 l200i 
iPuccetti et~all l2009j)_ using a m a ximum -likelihood rou- 
tine as detailed in ICivano et al.l (|2011ft . The Chandra 
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Fig. 1. — Rest-frame color U — V versus stellar mass for zCOS- 
MOS galaxies with 0.25 < z < 1.05 (small circles). Galax- 
ies in kinematic pairs are shown in black while all others are 
marked in grey. The AGNs are further marked with colored circles 
(filled=kinematic pair; open=all others). The vertical line marks 
our chosen mass limit. 
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Fig. 2. — X-ray luminosity (0.5-10 keV) versus redshift for AGN 
(large symbols) with those in galaxy pairs further marked with 
an open, red circle. The upper limit to the broad-band X-ray 
luminosity for each zCOSMOS galaxy is shown by the small grey 
dots. The solid horizontal line is the minimum X-ray luminosity 
that we enforce to measure the fraction of galaxies hosting AGN 
in two redshift intervals. 



observations cover the central 0.9 sq. deg and result in 
the detection of 1761 point sources above a flux level of 
5.7 x 1CT 16 erg cm" 2 s" 1 in the 0.5-10 keV band. The 
tiling scheme not only provides a fairly uniform sensitiv- 
ity but an optimal PSF across the field. We highlight 



that the resolving power of Chandra is important with 
respect to our study since galaxy pairs have been found 
to reside in reg i ons of heightened ga l axy overdensities 
(iLin et all [2010: Ide Ravel et all 1201 H iKampczvk et al.l 
1201 lh where extended X -ray emission, if present, may 
hamper the detection of point sources. The Chandra 
observations over the C-COSMOS region are sensitive to 
AGNs with Lo.g-io kcv > 2 x 10 42 ergs s~ 4 out to z ~ 1, a 
luminosity regime well-above that due to star formation 
thus indicative of a highly robust selection of AGNs. 

We identify 337 zCOSMOS galaxies with X-ray emis- 
sion (/0.5-10 kov > 1 x 10~ 15 erg cm -2 s -1 ) that pri- 



marily falls within the range 10 <Lq. 



5-10 keV 



< 10 4 



erg s _i (Figured]). The X-ray luminosity of each AGN 
is determined from their observed broad-band (0.5-10.0 
keV) flux assuming a power-law spectrum (photon index 
r = 2.0) and no correction for intrinsic absorption. In 
terms of Eddington ratio, we are sensitive to SMBHs hav- 
ing 10 -3 < LBoi/^Edd < 10 _1 given the stellar masses 
of their hosts, an assumption of their bulge-to-total mass 
ratio (B/T ~ 0.5), and the local Mbh ~ Mbuige relation. 
The majority of AGNs have optical emission dominated 
by their host galaxies and lack broad optical emission 
lines, both highlighting the importance of X-ray selection 
to account for the moderately-obscured population at 
these redshifts (jMainieri et al.ll2007t iBrusa et al.ll2010h . 

We compare the properties of our AGN hosts to the 
underlying galaxy population by selecting those with 



10 



42.0 



<Lo. 



5-10 keV 



< io 



43.7 



erg 



a range 



that 



likely results in a purely AGN-dominated sample, in 
terms of X-ray emission, and avoids the inclusion of 
more optically-luminous AGNs that can impact the de- 
rived rest-frame color s and stellar mass estimates (e.g., 
I Silverman et alJ l2008f l . In Figure [H we see that their 
host galaxies are massive (M* > 2.5 x 10 10 M Q ) and 
residual star formation may be present as evident by 
the slightly bluer colors (< U — V >= 1.6) of AGN 
hosts as compared to non-active galaxies (< U — V >= 
1.7) of similar stellar mass (e.g., [Silverman et alj|2009"bt 
ISchawinski et al.ll20 lOt fXue et al Jl20 ldh . This is in agree- 
ment with lKampczyk et alJ (|201lh who find that galaxies 
in zCOSMOS kinematic pairs have elevated star forma- 
tion rates, heightened post-starburst signatures and an 
increase in the numbers of galaxies with irregular mor- 
phologies. 

4. OBSERVED AGN FRACTION 

We first highlight some examples of AGNs associated 
with galaxies in close kinematic pairs (Figure [3]) that il- 
lustrate the diversity in their optical properties. In the 
top panel, a system of three galaxies at z ~ 0.659 is iden- 
tified, as shown in the HST/ACS F814W image (right; 
iKoekemoer et all [2007). to be kinematically linked. An 
X-ray source (left) is clearly associated with one of the 
galaxies undergoing an interaction, as is evident from 
the low-level optical emission linking it to a neighbor 
(towards the south-east) about 20 kpc away. In addi- 
tion to a point-like hard X-ray source, diffuse emission 
is present and more prominent in the soft energy band 
as expected (not shown). In Figure [3] (middle), two in- 
teracting spiral galaxies are separated by ~ 15 kpc with 
one of them hosting a luminous AGN while the other is 
a post-starburst galaxy based on spectral features (i.e., 
strong Balmer absorption series) evident in a zCOSMOS 
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Fig. 3. — AGNs in close kinematic pairs (left: Chandra; right: 
HST/ACS F814W). Top An AGN (CID=450) at z = 0.658 is part 
of a system of three galaxies within the central potential of an X- 
ray emitting galaxy group. The redshifts of the individual galaxies 
are labelled. Middle An AGN (CID=1711) is associated with a 
pair of interacting, spiral galaxies at z = 0.77. Bottom A system 
at z = 0.371 with a barred spiral hosting an AGN (CID=3083). In 
all panels, the optical positions of the zCOSMOS galaxies in pairs 
are marked with a small cross (red=AGN), and north is up while 
east is to the left. 



optical spectrum (not shown). A third example (bottom) 
clearly shows a barred spiral galaxy with an AGN and a 
companion 24 kpc away. 

Our objective is to measure the fraction of galaxies 
that host AGN as a function of whether they are asso- 
ciated with a kinematic pair or not in order to assess 
whether early encounters between galaxies induce such 
nuclear activity. We first perform an analysis based on 
"observed" galaxies in the zCOSMOS 20k catalog that 
are actually identified to be in the pair and the field 
samples. We note that the field sample includes galax- 
ies, in actual pairs, that have not been identified due to 
the spectroscopic incompleteness of zCOSMOS thus the 
observed AGN fraction in non-pair galaxies will likely 
be overestimated if an excess of AGNs in paired galax- 
ies is realized. In Section [5J we proceed a step further 
by recovering the intrinsic fraction of AGNs in galaxies 
in both pairs and the field, including reliable confidence 
intervals, since we have full knowledge of the selection 
function of the spectroscopic sample. This second step 
further enables us to compute the actual fraction of all 
AGN that are correlated with the presence of being in a 
kinematic pair. 

For all analysis, we specifically determine whether a 



zCOSMOS galaxy with M* > 2.5 x 10 10 M Q is identified 
with an X-ray detection by Chandra. We allow a galaxy 
to be associated with a kinematic pair if it has a neigh- 
bor satisfying the criteria, given above, on the projected 
separation and line-of-sight velocity difference. We fur- 
ther restrict the sample to pairs with a mass ratio less 
than 10 to 1 which allows for a companion to have a mass 
below our threshold. This mass ratio is chosen to isolate 
a sample for which strong gravitational interactions are 
capable of destabilizing gas that can subsequently fuel 
an AGN. 

The observed AGN fraction (/; equation [1]) is deter- 
mined by summing over the full sample of AGN (TV) 
with N e ff t i representing the number of effective galax- 
ies in which we could ha ve detected an AGN of a given 
X-ray luminosity L* x (see iSilverman et aT1l2009airbl ). 



N 



N, 



eff,i 



(1) 



This procedure enables us to account for the spatially- 
varying sensitivity limits of the Chandra observations (ef- 
fectively the source detection comp leteness as a function 
of X-ray flux shown in Figure 12 of iPuccetti et~aT1l2009l\ 
of th e COSMOS field (also see Figure 4 of lElvis et al.l 
2009). The effective number of galaxies typically falls 
below the total number of galaxies by ~ 15%. An ad- 
ditional weight (wi) denotes whether an X-ray source 
was observed as a random (wi = 1), or compulsory 
(wi = 0.70) target. This factor Wi is the ratio of the 
sampling rate of random targets (0.68) divided by that 
of the compulsory targets (0.97) within the central zCOS- 
MOS area. These sampling rates are simply the ratio of 
observed targets divided by the total number of targets 
in an input galaxy catalog. By down weighting the 'com- 
pulsory' targets, we effectively remove to first order the 
bias that some AGNs were targeted at a higher rate as 
compared to random galaxies. We also incorporate these 
weights (w) into N e ftj with a minor impact because the 
majority of galaxies (~ 97%) are random targets. Due 
to the fact that the Chandra observations were taken 
after the zCOSMOS program began, only 30% of these 
X-ray sources were compulsory targets, essentially the 
brighter XMM-Newton sources. We estimate the associ - 
ated la error (see equation 4 of iSilverman et al1 l2009b) 
based on a gaus sian approximation to a beta distribution 
(Cameron 2010, Andrae et al. in preparation). 

We then compare the fraction of galaxies, in pairs, that 
host an AGN to those, not in pairs, over the redshift 
range 0.25 < z < 1.05 (Table H). We define an AGN here 
as an X-ray point source with log Lo.5-10 kcV > 42.3 
(units of erg s _1 ), unless otherwise noted; the inclu- 
sion of a small number (13%) of more luminous AGNs 
(log Lx ~ 43.7 — 44.5), that may have less accurate stel- 
lar mass estimates, has no impact on our results. As 
motivated by the analysis given below, we limit the sam- 
ple of galaxies in pairs to those having a close neighbor 
that satisfies the criteria that dr < 75 kpc and dv < 500 
km s _1 . We find that over the entire redshift range the 
AGN fraction for the control sample is 5.0 ± 0.4% (not 
corrected for incompleteness with respect to pair identi- 
fication; see below), while the pair sample has a higher 
AGN fraction of 9.8 ± 2.2%. While the significance of 
the difference in the AGN fraction (A = f pa irs — f field) 
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Fig. 4. — Observed AGN fraction of galaxies in close kinematic 
pairs (dr < 75 kpc and dv < 500 km s" 1 ; filled black circles) as 
compared to galaxies with no neighbor within a projected separa- 
tion of 143 kpc and a velocity offset less than 500 km s — 1 (open 
circles) for two redshift intervals. The horizontal bars indicate the 
redshift range for each value while the vertical bars are the ltr 



is 2.2ct based on simple error propagation, we demon- 
strate below that the significance is actually higher when 
considering the sampling rates of the galaxy pairs (Sec- 
tion [5]). Even before taking this into consideration, an 
enhancement of AGN activity is evident for galaxies in 
pairs having close separations. 

In Figure |U we show the AGN fraction of galaxies, 
given this direct simple approach, in kinematic pairs and 
the control sample as measured in two redshift bins of 
equivalent width. We set the minimum X-ray luminosity 
to roughly the mean sensitivity limit for each redshift 
bin (see Table [T]) . We find that the fraction of galaxies in 
pairs hosting an AGN is higher in both redshift intervals, 
even when implementing a constant luminosity cut (i.e., 
log Lx > 42.6) over the full redshift range (not shown). 
This effectively demonstrates that our result based on 
the full redshift range is not due to a selection effect 
incurred by either the host galaxy mass limit, the mass 
ratio of the pairs, or an offset in the redshift distribution 
of AGNs and galaxies in pairs. The increase in the AGN 
fraction with redshift is a result of the overall evolution 
of the luminosity function of AGN. 

Our restriction on the projected separation (dr < 75 
kpc) of zCOSMOS pairs, in the above analysis, is due 
to the fact that we see an enhancement of AGN activ- 
ity specifically on smaller scales. To demonstrate this, 
we measure the AGN fraction of galaxies in pairs while 
separating the sample into two bins of projected sepa- 
ration (dr) and line-of-sight velocity difference (dv). In 
Figure Eh, we compare the results of this exercise with 
the AGN fraction of galaxies not in pairs that acts as our 
control sample. We find that the AGN fraction of galax- 
ies is highest (9.8%) for galaxies in pairs with dr < 75 kpc 
and a velocity difference of dv < 500 km s _1 . While the 
AGN fraction is elevated for AGNs in pair galaxies shown 



as a function of velocity difference (Figure [3)), the level 
of incidence is not arguably stronger in the bin with the 
smallest velocity offset. It does appear that a monoton- 
ically increasing AGN fraction (5.0%, 7.3%, 9.8%) with 
decreasing physical separation is evident when compar- 
ing the two dependent samples (dr < 75 kpc, dr < 143 
kpc) with the control (i.e., field) sample (first three lines 
in Table Q} . 

4.1. Mass ratios of galaxies in pairs 

We demonstrate that the bulk of AGN activity does oc- 
cur in zCOSMOS galaxies having a mass ratio of less than 
3:1 (i.e., 'major' interaction/merger). In Table[2] we give 
the AGN fraction of galaxies with selection in four dif- 
ferent mass ratio regimes (i.e., less than the given ratio) 
along with the size of the AGN sample. Our AGN selec- 
tion here includes a limit on the maximum luminosity of 
the AGN to avoid those that may have inaccurate stellar 
mass estimates. We find that 67% (12/18) of AGNs in 
galaxy pairs with dr < 75 kpc are associated with galax- 
ies having a mass ratio within this range. In Figure |6l we 
plot the AGN fraction as a function of mass ratios of the 
pairs and include all AGNs with log Lo.5-8 kev > 42.3. 
We remark that these measurements are not indepen- 
dent since they include all galaxies having a mass ratio 
less than the given value. There does appear to be a 
trend of increasing AGN fraction as the masses between 
the two galaxies are more similar. A larger sample is ob- 
viously required to justify such a claim with certainty. It 
is reassuring that for each mass ratio we find a consistent 
result that the AGN fraction increases in kinematic pair 
sample, as compared to non-pair galaxies (control sam- 
ple), with the highest levels of incidence occurring on the 
smaller scales of projected separation. We note that red- 
shift and mass dependent effects may be non-negligible 
for pair samples with higher mass ratios (> 10) since we 
do not implement a volume or mass limited criteria for 
both members of a galaxy pair. 

4.2. Further considerations 

We recognize that a number of selection effects must be 
considered given their potential impact on the aforemen- 
tioned results, such as (1) an artificially elevated level of 
incidence of AGNs in pairs due to differences in the red- 
shift success rate between galaxies and AGNs, or (2) 30% 
of the X-ray sources being spectroscopically-observed as 
compulsory targets. Any bias in the number of AGNs 
in our sample can artificially induce a higher fraction 
of pairs with AGNs if not considered carefully since the 
probability of having a pair is the product of the sam- 
pling rates of the individual galaxies. To address the first 
issue, we measure the success rate of r edshift identifica- 
tion with confidence flag / > 1.5 (see iLillv et al.l [20091) 
between randomly-targeted galaxies (89.2%) and AGNs 
(91.9%) in the zCOSMOS sample and find very simi- 
lar rates. To address the second issue, we carried out 
a similar analysis, as done above, to measure the AGN 
fraction but with a highly conservative selection. We lim- 
ited the sample to (1) only galaxies and AGN that were 
observed as part of the random sample (see Section [5]). 
In addition, we constructed a highly conservative sam- 
ple that further imposed the following restrictions on the 
random sample: (2) only AGNs with a higher luminos- 
ity cut (log Lx > 42.6) over all redshifts to avoid those 
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Fig. 5. — Observed fraction of galaxies hosting an AGN shown as a function of both the projected physical separation (a) and line-of-sight 
velocity difference (6). The right-most data point in both panels gives the AGN fraction of galaxies not in pairs with the horizontal lines 
showing the lcr error and extended along the abscissa for visual comparison. 



TABLE 1 

Sample statistics-observed kinematic pair analysis 3. 



Type 


Redshift 


dr 


dv 




#of 


#of 


AGN 






(kpc) 


(km s" 1 ) 




AGN 


galaxies 


fraction (%) 


Pairs 


0.25-1.05 


< 75 


< 500 


>2 


19 


223 


9.75 ±2.15 




0.25-1.05 


< 143 


< 500 


> 2 


38 


562 


7.30 ±1.17 


Control 


0.25-1.05 


> 143 


> 500 


> 2 


116 


2726 


5.03 ±0.45 


Pairs 


0.25-0.65 


< 75 


< 500 


> 2 


7 


102 


6.75 ±2.56 




0.65-1.05 


< 75 


< 500 


> 4 


14 


138 


11.44 ±2.96 


Control 


0.25-0.65 


> 143 


> 500 


> 2 


23 


907 


2.56 ±0.54 




0.65-1.05 


> 143 


> 500 


> 4 


87 


1819 


5.24 ±0.56 


Pairs-random 


0.25-1.05 


< 75 


< 500 


> 2 


14 


214 


8.25 ±2.05 










4-50 


13 


214 


7.10 ±1.88 


Control-random 


0.25-1.05 


> 143 


> 500 


> 2 


87 


2671 


4.28 ±0.43 










4-50 


63 


2671 


2.70 ±0.33 



a Not corrected for incomplete identification of galaxy pairs 
b units of 10 42 erg s _1 (0.5-10 keV) 



the case where the X-ray luminosity is limited to the 
lower end of the distribution (42.6 < log Lx < 43.7) 
that may indicate that such encounters are more likely 
to result in AGN of moderate luminosity. 

It may also be the case that the observed trends are 
driven by the dependence of AGN activity on stellar 
mass. Therefore, we performed a series of K-S tests to de- 
termine whether the stellar-mass distribution of the pairs 
and the control sample differ in any of the bins (i.e., red- 
shift, dr, dv) implemented in this study. We find that 
there is no noticeable difference that can be responsible 
for the results presented herein. 

5. A MONTE-CARLO INFERENCE OF THE 
INTRINSIC AGN FRACTION 

As a second step, a complementary analysis, using a 
Monte Carlo technique as fully outlined in the Appendix, 
is implemented to recover the intrinsic fraction of AGNs 
in different subsets of galaxies such as those in galaxy 
pairs or within a field sample. This procedure allows us 
to properly account for (1) the selection effects outlined 



TABLE 2 
Fraction of galaxies 3 in pairs of 
different mass ratios hosting agn b 



Mass ratio c 


dr < 75 kpc 


dr < 143 kpc 




% (# AGN) 


% (# AGN) 


15:1 


9.3 ±2.1 (19) 


6.6 ± 1.1 (35) 


10:1 


9.8 ± 2.2 (18) 


6.7 ± 1.1 (33) 


5:1 


11.1 ±2.6 (17) 


7.4 ± 1.3 (29) 


3:1 


11.7±3.2 (12) 


7.6 ± 1.6 (21) 



a Not corrected for incomplete identification 
of galaxy pairs 

b log I/O. 5-10 keV > 42.3 

c Mass ratio less than the given range. 



falling near the flux limit thus avoiding a large contri- 
bution from a few objects, (3) as done above, a limit on 
the maximum luminosity of the AGN. Even with these 
limited samples, we see an enhancement of AGNs on the 
same scale and with similar confidence (see Tabled]). We 
highlight that a stronger enhancement is seen (2.6 x) for 
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Fig. 6. — Observed AGN fraction of galaxies in close kinematic 
pairs (filled circles: dr < 75 kpc and dv < 500 km s — 1 ; open circles: 
dr < 143 kpc and dv < 500 km s" 1 ) as a function of the mass ratio 
of the two galaxies. The horizontal bars give the non-pair fraction 
with lcr errors. 



above, (2) the incompleteness with respect to pair iden- 
tification (not yet addressed), and (3) determine confi- 
dence intervals for all quantities without an assumption 
such as Gaussianity. Accurate knowledge of the fraction 
of galaxies in a non-pair sample enables us to estimate 
the impact of mechanisms, not associated with external 
interactions (i.e., secular processes), responsible for AGN 
activity and provide an accurate determination of the 
boost factor of AGNs in pairs relative to the field sample, 
i.e. the excess of AGN supposedly induced by being in a 
kinematic pair, after taking out the field "background". 

The results of this exercise are shown in Figure [7] and 
given in Table [3] (including 68% and 99% confidence in- 
tervals). We find a consistent result with the previous 
analysis for which the AGN fraction of galaxies in pairs 
with dr < 75 kpc is 9.7% (see panel c). It is now ev- 
ident that the AGN fraction in the field is significantly 
reduced to 3.8%. Based on these values given in Table |3] 
and their confidence intervals, we now find an intrinsic 
enhancement of AGNs in pairs by a factor of 2.6 with 
a higher significance (3.2c) than previously determined. 
Although, an accurate assessment of the significance is 
non-trivial given the anti-correlation between /p a j rs and 

Afield (Fig. [TJi, b). Also consistent with the previous anal- 
ysis, we find an increase in the AGN fraction in pairs, 
spanning the full range of separation (dr < 143 kpc; 
Fig.[7]i), as compared to the field estimate, although not 
as strong an enhancement as seen on smaller scales. 

5.1. Fraction of the AGN population attributed to 
galaxy interactions 

We can now use our measurement of the intrinsic frac- 
tion of galaxies hosting AGN, either associated with 
kinematic pairs or not, to determine the relative im- 
portance of galaxy interactions, as compared to other 
factors yet to be firmly established. These other pro- 



cesses are likely to be internal to their host galaxy since 
we have determined, using the same zCOSMOS galax- 
ies, that AGN activity is also depende nt on the stellar 
mass a nd star formation rate of its host ()Silverman et al.l 
2009b) and not stro ngly dependent on the large-scale 
environment (~Mpc; Silverman et al.1 l200 9a) . Further, 
iCisternas et al.l ()2011 ) demonstrate that the majority 
(80%) of AGN hosts in COSMOS are not undergoing ma- 
jor mergers. We highlight that the enhancement of AGNs 
in pairs is really a local effect (< 143 kpc); the fact that 
zCOSMOS pairs are more common in regions of height- 
ened large-scale galaxy overdensity (jde Ravel et al.ll20lTl : 
iKampczyk et al.l 1201 If ) does not mean that structures 
such as galaxy groups are more conducive for AGN ac- 
tivity (e.g., through the accretion of cold gas) beyond in- 
creasing the galaxy merger rate. We specifically find that 
AGNs are equally likely to reside in optically-selected 
grou ps and the 'field' (non-g roup) population (see figure 
11 of lSilverman et al.ll2009af ) . 

To address this issue, we measure an excess fraction 
of galaxies in pairs as compared to those in the 'field'. 
This ratio (f x ) is the number of AGN triggered by in- 
teractions (n x ) relative the the total AGN population 
(n 



pairs 



l ficld' Equation [5]) . The number of AGN in 

pairs (n x ) can be expressed with n^g (bg=background), 

the number of AGNs in the pair sample not triggered by 
interactions. 



pairs 



fx 



n x 



"x 



"pairs + "field 



(2) 



l bg 



n x 



/field x ^pairs ( 3 ) 



These equations can be combined and arranged to give 
the following expression. 



fx 



pairs 



/field 



/, 



pairs 



/field 



'field 
pairs 



(4) 



We can then evaluate this expression through the re- 
sults of our Monte Carlo approach to measure the intrin- 
sic AGN fraction in the field and in pairs (see Table ■ 
In the lower two panels of Figure[71 we plot the likelihood 
distribution of f x for two pairs samples (left: dr < 75 
kpc; right: dr < 143 kpc). We find that interactions be- 
tween galaxies with dr < 143 kpc account for 17.8^74% 
of AGN activity in our sample. 



6. CONCLUDING REMARKS 

We have performed a simple test of the scenario for 
which galaxy interactions trigger moderate-luminosity 
AGNs. To do so, we have utilized a sample of kine- 
matic pairs identified from the zCOSMOS 20k 'bright' 
catalog and Chandra observations that indicate those 
harboring AGNs. The X-ray selection of AGN enables us 
to include those that may be cloaked in star formation 
that can hamper optical selection, especially since ele- 
vated SFRs are a common property of close galaxy pairs 
(|Ellison et alJl20"08t ILin et al.ll2010l : ILiu et alJl2lHll) . 

Based on the multi-wavelength observations in COS- 
MOS, we have found that galaxies (0.25 < z < 1.05; 



Silverman et al. 



dr <75 kpc 



dr<U3 kpc 




0.02 0.04 0.06 0.08 0.10 0.12 0.14 

AGN fractions 



0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 

AGN fractions 




0.00 0.02 0.04 0.06 0.08 0.10 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

^ f vairs f field ^ f pairs f field 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 

P J pairs f field p f pairs J 'field 

jx = jzr jx ~ 



J pairs • J field 



J pairs ' J field jy 



Fig. 7. — Intrinsic fraction of galaxies hosting AGN based on a full Monte Carlo analysis. The rows are as follows: (1) AGN fraction 
in the field versus pair sample for the individual realizations, (2) probability distribution for the AGN fraction in the field (blue) and pair 
sample (red), (3) difference in the AGN fraction between pairs and the field, (4) excess fraction as defined in the text. The left column 
refers to pairs with dr < 75 kpc while the right column shows the equivalent for pairs with dr < 143 kpc. 



TABLE 3 

Intrinsic kinematic pair analysis (Monte Carlo approach) 



Type Redshift dr dv L x a AGN fraction (%) Excess fraction (f x ) 

(kpc) (km s- 1 ) (1%, 16%, median, 84%, 99%) (1%, 16%, median, 84%, 99%) 



Pairs 0.25-1.05 < 75 < 500 > 2 6.3, 8.0. 9.7, 12.0, 15.6 6.1, 11.9, 17.6, 24.1, 34.0 

0.25-1.05 < 143 < 500 > 2 4.8, 5.8, 6.7, 7.7, 8.9 2.2, 10.4, 17.8, 26.2, 35.1 

Field 0.25-1.05 > 143 > 500 > 2 3.0, 3.4, 3.8, 4.1, 4.5 — 



a units of 10 42 erg s" 1 (0.5-10 keV) 

M* > 2.5 x 10 10 M Q ) in kinematic pairs, with physi- 
cal separations less than 75 kpc and line-of-sight veloc- 
ity difference less than 500 km s _1 , exhibit statistically 
significant (3.2er), intrinsically heightened levels of AGN 
activity (9.7j^ - 7%) relative to galaxies not within these 
associations (3.8^ '|%) by a factor of 2.6. Interestingly, 
this factor is remarkably close to the boost in star forma- 
tion seen in zCOSMOS kinematic pairs on small scales 
(jKampczvk et all 120 111) thus possibly indicating a com- 



mon causal connection between galaxy interactions, star 
formation and AGN activity. 

Our observed enhancement of AGNs on these phys- 
ical scales (dr < 75 kpc) is in broad ag reement with 
merger-driven models (jHopkins et al.ll2~008|) of black hole 
growth. From these models, the first signs of an ob- 
scured AGN should emerge roughly 0.5 Gyr after first 
passage of two interacting galaxies and of similar physi- 
cal separation. Based on a comparison to the Millenium 
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Fig. 8.— HST/ACS images of AGNs that fall within our field 
sample and identified to be in late stages of a galaxy merger based 
on their optical morphology of the host galaxy: (a) CID799, (b) 
CID298, (c) CID401, (d) CID496 

mock catalogs of the COSMOS field, iKampczvk et al.l 
(|2011l ) demonstrate that the majority of zCOSMOS kine- 
matic pairs, with such small separations, are likely to 
merge within a Hubble time. Keep in mind that our 
results pertain mainly to moderate-luminosity AGNs 
(Lx ~ 10 43 erg s" 1 ), and in no way are we claiming that 
kinematic pairs are conducive environments for the more 
luminous quas ars. Although, it has been shown that 
SDSS quasars (Scrbcr ct al. 2006) are more strongly clus- 
tered on scales less than 100 kpc. 

There has been much debate as to whether or not 
galaxy interactions including major mergers are the sin- 
gle mechanism predominantly driving gas to the nuclear 
region in galaxies thus inducing accretion onto the central 
SMBH. It is now apparent, based on ma ny studies (e.g., 
iGrogin et all 120051 : iGabor et al.l I2009L ICisternas etaLl 



1201 If) including the one presented here, that a number of 
physical processes, yet to be fully understood, are impor- 
tant. Therefore, we can now attempt to determine the 
relative contribution of various mechanisms responsible 
for black hole growth. Based on our results, we start by 
estimating the fraction of the moderate-luminosity AGN 
population that can be attributed solely to galaxy inter- 
actions; we find that ~ 18% of such AGNs are the result 
of interactions. This leaves open many questions such as 
the following: What physical mechanisms are responsible 
for the remaining ~ 80% of AGNs? We do recognize that 
late-stage mergers are not accounted for in our sample. 
It is clear that these do exist within the field population 
(Figure with some exhibiting double optical nuclei, 
although it is unlikely tha t they account for a m ajority 
the AGN population since (|Cisternas et alJfeOllh showed 
that at z < 1 the fraction of moderate-luminosity AGN 
that are due to the late stages of interaction has an up- 
per limit of ~ 20%. Whether the final coalescence of two 
massive galaxies results in a luminous quasar remains 
to be determined and if they emerge at a later times in 
merger sequence. Answers to these questions may require 
further Chandra observations of the full zCOSMOS area 
or await larger area spectroscopic surveys having a high 
degree of completeness. 
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APPENDIX 

MONTE CARLO SIMULATION OF INTRINSIC AGN FRACTIONS 

We construct a method to infer the intrinsic AGN fraction of galaxies either in pairs or the field based on the 
observed galaxy sample with full consideration of the dominant selection effects. The analysis is complicated by two 
special characteristics of our data. First, galaxies have a detection probability of w « 0.60, irrespective of whether 
they host an AGN or not, due to the spectroscopic completeness of the zCOSMOS 'bright' survey (see iLillv et al.1 
120071 [2009h . Here, we define spectroscopic completeness as the ratio of galaxies (jab < 22.5) with secure spectroscopic 
redshifts (confidence flag > 1.5) to all galaxies above this magnitude in our photometric catalog. As a consequence, 
galaxies that actually reside in pairs, but were not identified as such (since the partner was not observed), will be 
present in the field sample. Second, there is a subsample of Nqs compulsory targets, which are known a priori likely to 
contain an AGN based on their X-ray or radio emission, although the former is of most concern for our study. These 
compulsory targets have detection probability of 1, in contrast to the other "random" targets. These complications 
have to be addressed in order to obtain an unbiased estimate of the intrinsic AGN fractions. While we do account for 
the differential weighting (compulsory versus random) between targets in Section IH the fact that galaxies in actual 
pairs can be present in the field sample is not addressed in that analysis. 

To fully account for these complexities with the data, we set up a probabilistic model based on a decision tree that 
accounts for all possible outcomes of an observation. Given an intrinsic total number of galaxies that could have 
been observed down to iab = 22.5 within the joint zCOSMOS 'bright' and Chandra survey area, the decision tree 
considers the following possibilities: 

• Is the target random or compulsory? 

• Has the target been observed? 

• Is the observed target an AGN? 
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• Does the observed target reside in a pair? 

• If the target resides in a pair, has the partner been observed? 

The decision tree has free model parameters that have to be inferred from the data. Among these model parameters 
are the desired intrinsic AGN fractions in the held // (dr > 143 kpc and dv > 500 km s -1 ), in close pairs fa (dr < 75 
kpc and dv < 500 km s" 1 ) and in loose pairs /143 (75 < dr < 143 kpc and dv < 500 km s" 1 ). Other parameters are 
the intrinsic total number of galaxies A^^, and the probability of a galaxy or AGN to reside in a close pair p 7 $ or in 
loose pairs P143. The decision tree then predicts the observed quantities. Introducing the number of random targets 
A^ran = ^ot ~ ^ cs and * ne probability of a target not to reside in a pair but in the field pf = 1 — P75 — P143, the 
predicted number of observed galaxies in the field is 

ArP red = (N CS + N T& nw) (p f + (p 75 +p 143 )^™£(l - w ) j (Al) 

V ^tot / 

and the predicted number of observed AGN in the field is 

n P red = N cs ( p f + (p 75 +Pi43) ; ^j(l - w) J + Aranw (pfff + (^75/75 + Pi43/i43) ; ^^jr(l - w) J . (A2) 
V ^tot / V ^tot / 

Similarly, the predicted number of observed galaxies in close or loose pairs is 



pred / at , \ ( N cs . N Tan \ 

V v tot Jv tot/ 



<5/143 = P75/143 (^CS + wN raTl ) ( ^ + W-££ ) (A3) 



and the predicted number of observed AGN in close or loose pairs is 

pred 



75/143 



P75/i43(^cs + w ■ / 7 5/i43^ran) (A^cs + A ran w) . (A4) 



The free parameters of the decision tree, A^ n jj , p 7 $ , P143 , ff, fa, and /143, then have to be estimated in order to match 
the observed numbers. 

In fact, for given observed number of galaxies and AGN in field and pairs, the set of equations given by Eqs. (|Aip . 
(|A2[> . (|A3|) . and (|A4I) can be solved numerically. These numbers provide an initial guess for the free parameters of 
the decision tree. While the central-limit theorem suggests that the numerical solution is a good approximation to 
the most likely values for our parameters, we are also interested in evaluating their uncertainties. Accurate error 
estimation is crucial to determine whether an excess of AGN activity in galaxy pairs is realized. In order to obtain 
confidence intervals on our free parameters, we need to introduce a likelihood function to our probabilistic model as 
given by 

r( AT obs Ajobs i\robs „obs „obs ^obsiArint _ f f f \ _ 

L(A 7 g ,A 143 ,A f ,n 75 ,ra 143 ,rc f |A tot ,P75,Pl43, //, /75, /143) - 

= C(N^ S \N^l,p 75 ,p 143 , ff, fa, fa 3 ) ■ C(rff | N™\ , p 75 , Pl43, //, fa, /l4 3 ) 

• C(N$§\N™l,p 75 ,p li3 , ff,fa, A43) • £(n^\N$lp 7n ,p U3 , ff,fa, fas) 

■ £(N? hs \N$l,p 75 ,p 143 , f f , fa,f 143 ) ■ C(nf s \Nlf t ,p 75 ,p 143 Jf,fa, f 143 ) , (A5) 

which can be factorized because the measurement processes of N^ s , ^^43) N^° s , n^g s , , an d n 1 ^ 08 are statistically 
independent of each other. Of course, the underlying astrophysics imposes relations of these observables as described 
by our probabilistic model. However, these relations do not lead to correlated measurements. The model predictions 
of Eqs. (|A1IA2IA3IA4[) are always the sums of two binomial variates s = b\ + b 2 , whose likelihood is then given by the 
convolution of the two binomial likelihoods p(bi\Ni, fa 

/>min(s,7Vi) 

C(s\N 1 ,N 2 ,faf 2 )= p{b\N x ,h)p{s-b\N 2 J 2 )db. (A6) 

As fi ^ f 2 in our case, this convolution integral can only be solved numerically not analytically, which makes the 
parameter estimation computationally expensive. For instances, from Eq. (|A3[) . we can identify the parameters 



/ N cs , A ran \ / N C8 N mn \ ( , 7 , 

h = P ^ 735t + W ^hrt and h = P ™/^ W 735t + W ^hrt ■ < A7 ) 
V v tot Jv tot/ V v tot Jv tot/ 



This enables us to evaluate the likelihood function Eq. (IA5|) of the data for given values of the model parameters 
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Having identified the likelihood function of the observed data, we can now estimate the model parameters by 
maximizing this function. Again, since we are particularly interested in assessing differences between the intrinsic 
AGN fractions, reliable error estimation is of vital importance. Maximum-likelihood parameter estimation is favorable 
in this case because it provides minimal uncertainties on the results (e.g., Barlow 1998). Other routines for parameter 
estimation are applicable but inevitably produce estimates with larger errors. One option is to draw on the central-limit 
theorem and approximate the likelihood function by a Gaussian at its maximum. Error estimates can then be derived 
by Fisher analysis (Heavens 2009). However, we do not want to rely on any such approximations, given the number 
statistics, especially due to the low AGN fraction of our sample. Sin ce our probabilistic model contains several free 
parameters and the evaluation of the convolution integral in Eq. (IA6[) is expensive, we employ a Metropolis-Hastings 
algorithm since a grid search is not an option. Using such a Markov-Chain Monte-Carlo algorithm has the additional 
advantage of easily providing marginal estimates /143, /f, and their differences, which arc independent of the other 
model parameters. This marginalisation is necessary, because otherwise we only obtain conditional error estimates 
that underestimate the actual errors and would lead us to report overly decisive differences in AGN fractions. 
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